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Simulations of a supersonic recessed-cavity flow are performed 
using a hybrid large-eddy/Reynolds-averaged simulation 
approach utilizing an inflow turbulence recycling procedure 
and hybridized inviscid flux scheme. Calorically perfect air 
enters a three-dimensional domain at a free stream Mach 
number of 2.92. Simulations are performed to assess grid 
sensitivity of the solution, efficacy of the turbulence recycling, 
and the effect of the shock sensor used with the hybridized 
inviscid flux scheme. Analysis of the turbulent boundary layer 
upstream of the rearward-facing step for each case indicates 
excellent agreement with theoretical predictions. Mean 
velocity and pressure results are compared to Reynolds- 
averaged simulations and experimental data for each case and 
indicate good agreement on the finest grid. Simulations are 
repeated on a coarsened grid, and results indicate strong grid 
density sensitivity. Simulations are performed with and 
without inflow turbulence recycling on the coarse grid to 
isolate the effect of the recycling procedure, which is 
demonstrably critical to capturing the relevant shear layer 
dynamics. Shock sensor formulations of Ducros and Larsson 
are found to predict mean flow statistics equally well. 
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1. Introduction 

In order to address some of the fundamental limitations of 
Reynolds-averaged simulation (RAS) approaches while also 
exploiting the strengths of large-eddy simulation (LES) 
approaches, researchers frequently use a blended LES/RAS 
approach. This is especially true for high Reynolds number 
supersonic wall-bounded flows, which exhibit inherent 
unsteadiness, turbulent shock-boundary layer interactions, 
and other complex flow phenomena. Whereas the RAS 
equations require modeling of turbulence at all length 
scales, the LES equations only require modeling of turbulent 
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motions at length scales below a characteristic filter width; 
the largest turbulent motions are resolved naturally through 
the solution of the filtered governing equations. Since these 
larger turbulent structures are responsible for the majority of 
mass, momentum, and energy transport, LES simulations 
are well-suited for high Reynolds number flows. 
Unfortunately, the grid densities required for LES 
simulations of wall-bounded flows often render the 
approach computationally infeasible for problems of 
practical importance. 

In this study, a hybrid LES/RAS approach is used in 
conjunction with a turbulence recycling procedure to 
simulate a supersonic flow over a rearward facing step and 
subsequent reattachment along an inclined wall [1-3]. The 
recycling process sustains coherent turbulent structure 
within the inflow boundary layer without having to simulate 
its development. Simulation results are compared to 
experimental mean flow data within the approach turbulent 
boundary layer, through the shear layer, and along the 
inclined wall. Assessments of the LES/RAS simulations are 
made using steady- state RAS results and experimental data. 

The simulations are designed to assess several recent 
developments to the Viscous Upwind aLgorithm for 
Complex flow ANalysis (VULCAN) computational fluid 
dynamics (CLD) code, including: reduced-dissipation 

numerics, a hybridized inviscid flux scheme, and inflow 
turbulence recycling [4,5]. Specifically, the current study: 
assesses the applicability of the hybrid LES/RAS method 
for cavity flows, isolates the effect of inflow turbulence 
recycling on the flowfield statistics, determines any grid 
sensitivity, and appraises the choice of shock sensor for use 
with the hybridized inviscid flux scheme. These objectives 
are achieved by performing seven simulations, which are 
parameterized by solution method, inflow treatment, grid 
size, and shock sensor and are summarized in Table 1. 
Hereafter, each simulation will be referenced using the case 
ID stated in Table 1. Lor example, references will be made 
using “case R35D”, rather than “case one”, where R, 35, 
and D indicate recycling, 35 million grid cells, and Ducros 
sensor, respectively. 

In the following section, the physical flow and wind-tunnel 
experiment are described. In section three, the numerical 
formulation is presented. After describing the grids and 
computational execution in section four, the results and 
accompanying discussion are presented in section five. 
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Conclusions and summarizing comments are included in the 
final section. 


Table 1 — Simulation cases parameterized by inflow 
treatment, grid size, and shock sensor. 


Case ID 

Solution 

Method 

Inflow 

Treatment 

Grid 

Cells 

Shock 

Sensor 

R35D 

LES/RAS 

Recycling 

35e6 

Ducros 

R35L 

LES/RAS 

Recycling 

35e6 

Lars son 

R8D 

LES/RAS 

Recycling 

8.8e6 

Ducros 

R8DW 

LES/RAS 

Recycling 

35e6 

Ducros 

R8L 

LES/RAS 

Recycling 

8.8e6 

Lars son 

NR8D 

LES/RAS 

No 

Recycling 

8.8e6 

Ducros 

NR8RAS 

RAS 

No 

Recycling 

8.8e6 

None 


2. Physical Flow 

The flow chosen for the current investigation was first 
studied by Settles, et. al. in the early 1980s [1-3]. The 
original experiment was performed in the Princeton 
University 20cm by 20cm High Reynolds Number 
Supersonic Wind Tunnel and consisted of Mach 2.92 air 
entering the tunnel test section at a stagnation temperature 
of 25 8 K, a stagnation pressure of 0.69MPa, and a free 
stream unit Reynolds number of 6.7x1 0 7 /m. The test 
article, presented in Fig. 1, was designed such that a 
turbulent boundary layer developed along a flat plate before 
separating over a sharp corner. The free shear layer formed 
over the subsequent cavity and reattached along a wall of 
length 18cm inclined at 20°. The position of the inclined 
wall was adjusted before taking measurements until no 
change in pressure or flow direction occurred due to the 
separating boundary layer, and the experimenters observed 
that the shear layer exhibited significant sensitivity to the 
inclined wall positioning. The cavity depth was 2.54cm, 
and the cavity assembly was unaffected by the tunnel 
sidewall boundary layers due to its positioning away from 
the walls. In order to promote two-dimensionality, 
aerodynamic fences lined the edges of the inclined wall. 
Mean flow measurements were made using a traversing set 
of pitot and static pressure probes, and these measurements, 
in conjunction with an assumption of constant total 
temperature, enabled the researchers to calculate Mach 
number, velocity, density, and other properties. A spark 
shadowgraph enabled observation of major flow features 
(not shown in this paper). 

The original experiment was designed to illustrate the 
complex fluid phenomena comprising the reattachment of a 
free shear layer in a supersonic flow. The effect of 
turbulence on the flowfield is prominent, as documented in 
the original review of the experiment and in several 
subsequent numerical investigations [2,3]. Major flow 
features observed in the experiment are illustrated in Fig. 2. 


After the approach boundary layer detaches at the step 
corner, a free shear layer forms over the cavity. As the 
shear layer reattaches at the inclined wall, fluid is entrained 
into the cavity, thereby forming a recirculation zone. The 
inherent unsteadiness of the shear layer alters the location of 
reattachment and subsequently affects the fluid entrainment. 
Together, the shear layer unsteadiness, altered attachment 
point, and recirculation zone form a feedback loop that give 
rise to a low-frequency flapping motion of the shear layer. 
Furthermore, as the flow reattaches, an oblique compression 
shock forms off the inclined wall. 



Figure 1 — Experimental setup with one aerodynamic 
fence removed and the current computational domain 
footprint signified by dashed lines. 



Figure 2 — Centerline profile of the experimental 
apparatus with major flow features highlighted. 

3. Numerical Formulation 

All simulations are performed using the VULCAN CFD 
code, which solves the Navier-Stokes equations using a cell- 
centered finite- volume approach. Hybridization of the LES 
and RAS equations is achieved through the blending of the 
respective closure models. Accordingly, the user has 
significant flexibility in choosing both the RAS turbulence 
model and the LES subgrid-scale (SGS) model. 
Development of the hybrid LES/RAS approach is provided 
by Baurle, et. al. [6], and details of the current 
implementation are described by Baurle, et. al. [7]. 

For each of the LES/RAS simulations, the Menter-BSL 
turbulence model is used to close the RAS equations [8]. 
This model requires the solution of two additional 
equations — a transport equation for turbulent kinetic energy 
and a transport equation for turbulence frequency. To 
model the SGS turbulence for use in LES regions, the one- 
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equation model of Yoshizawa and Horiuti is used and 
requires the solution of a transport equation for the SGS 
turbulence kinetic energy (TKE) [9]. 

A spatial- and flow-dependent blending function is used to 
blend the LES and RAS models throughout the 
computational domain. As designed, the current blending 
function transitions the RAS formulation to the LES 
formulation in the outer region of the boundary layer, near 
the upper edge of the log layer. As previous studies have 
indicated [10], targeting the transition for this region at the 
edge of the log layer within the turbulent boundary layer 
yields the most resolved turbulence without adversely 
affecting the universality of the boundary layer. The 
blending function used in the current study is presented by 
Edwards, et. al. [11]- 

In order to reduce numerical dissipation in smooth flow 
regions, while maintaining the dissipation necessary to 
capture discontinuities elsewhere, a hybridized flux scheme 
is used to construct interface inviscid fluxes. This 
hybridization is implemented as the linear combination of a 
low-dissipation, symmetrically-reconstructed inviscid flux 
and a conventional upwind-biased flux. The low-dissipation 
flux formulation uses left and right states constructed to 
higher-order accuracy at the cell interfaces using a four- 
point symmetric interpolation on the primitive variables 
[12]. A dissipative flux is then computed by first using the 
piecewise-parabolic method of Collela and Woodward to 
construct primitive variables at either side of a cell interface 
and subsequently using the Harten, Lax, van Leer, and 
Contact (HLLC) scheme to compute the dissipative flux in 
conjunction with a minmod-type limiting procedure [10,13- 
15]. The final hybridized inviscid interface flux is 
computed by blending the low-dissipative and dissipative 
fluxes using a shock sensor weighting parameter. The two 
shock sensor formulations of Ducros [16] and Larsson [17] 
are used in this study. Each of the shock sensors are 
designed to dynamically adjust the inviscid flux calculation. 
In regions of viscous, smooth flow, the shock sensor 
emphasizes the low-dissipation formulation more heavily, 
whereas in regions of large inviscid gradients, the 
dissipative flux formulation is more heavily weighted. The 
lower bound for each of these sensors was modified to 
provide a small level of background dissipation for the 
current work. Rather than allowing the sensor to range from 
0.0 to 1.0, each shock sensor is bounded by 0.2 and 1.0, 
thereby preventing numerical instabilities arising within the 
shear layer reattachment region. A detailed development of 
this hybridized method is presented by White, et. al. [5]. 

In order to sustain coherent structure within the turbulent 
inflow boundary layer, turbulent fluctuations and relevant 
turbulence properties are recycled from a plane one cavity- 
depth upstream of the step. Fluctuations are extracted from 
the donor plane by subtracting time- and span-averaged 
mean flow data from the instantaneous flow data. These 
fluctuations are scaled according to boundary-layer 
similarity laws before being superimposed to a mean flow 


profile. Initial fluctuations are provided by a previously- 
simulated fully-resolved (via LES) turbulent boundary layer 
on a flat plate, which provides the turbulence necessary for 
efficient recycling without necessitating the natural 
development of such structures. Details of the recycling 
procedure are described by Choi, et. al. [10], Xiao, et. al. 
[18], and Fan, et. al. [19]. 

4. Computational Execution 

Simulations are performed on structured three-dimensional 
fine and coarse grids of 35 and 8.8 million cells, 
respectively. The coarse grid was constructed by removing 
every other grid point in each of the streamwise and 
spanwise directions of the fine grid. The grid generated for 
this study is composed of four primary blocks, as shown in 
Fig. 3. The two upper blocks are connected to the lower 
blocks via non-CO continuous patches. Similarly, the lower 
left block is attached to the lower right block by a non-CO 
continuous patch where only the spanwise coordinates are 
mismatched (the downstream spacing is twice that of the 
upstream spacing). On the fine grid, streamwise-trans verse- 
span wise dimensions of the blocks adjacent to the wall 
directly upstream and downstream of the step are 
525x121x301 and 453x241x151, respectively; dimensions 
of the upper left and right blocks are 29x45x15 and 
153x45x15, respectively. The patches connecting the lower 
and upper blocks are positioned far enough from the domain 
of interest to ensure that resolved boundary layer and free 
shear layer eddies are not convected into them. 



Figure 3 — Three-dimensional domain coarsened twice in 
the streamwise and transverse directions and once in the 
spanwise dimension for visual clarity. 

The physical dimensions of the computational domain are 
28.6cm, 22.5cm, and 3.81cm in each of the streamwise, 
transverse, and spanwise directions, respectively. The flat 
plate extended four step heights upstream of the step, which 
is a factor of two further upstream than in previous studies 
to allow for recycling of turbulent fluctuations. The grid 
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was clustered in the wall-normal direction to provide a 
nominal y+ of 1.0 [5,20]. The stream wise and span wise 
grid spacing was set to l/20th of the boundary layer 
thickness. Downstream of the step, the grid clustering 
required by the approach flow boundary layer was relaxed 
toward a more isotropic grid for the resulting shear layer. 
The average streamwise, transverse, and spanwise spacing 
in this region was 2.0%, 0.2%, and 1.0%, respectively, of 
the cavity depth (2.54cm). The cell aspect ratio for the LES 
regions was no greater than approximately 20.0. The wide 
grid used for case R8DW extended 4.5 cavity depths in the 
spanwise direction and maintained the same grid spacing 
and clustering as the coarse grid. The width of the 
computational domain was increased by a factor of three to 
investigate potential correlations between two-point 
statistics that may be present in the original domain size. 
The three-dimensional fine and coarse grids were split into 
1980 sub-blocks to allow runtime parallelization using 
Message Passing Interface (MPI) on NASA Advanced 
Supercomputing resources. Each simulation performed on 
the fine grid required approximately 360 hours of wall clock 
time on 512 processors. The same simulations performed 
on the coarse grid required about 90 hours of wall clock 
time on 512 processors. 

Inflow profiles and flowfield initializations for the 
LES/RAS cases were obtained using the NR8RAS case, in 
which the Menter k-co turbulence model was used with a 
compressibility correction [8,21]. Further, the free stream 
flow was allowed to develop until the momentum thickness 
one cavity-depth upstream of the step matched experimental 
measurements. A significant spread in solutions was 
obtained using various turbulence models for the 
preliminary RAS simulation, and the RAS results used in 
this study represent those in best agreement with 
experimental data. Nominal inflow conditions included: 
static pressure of 21.2kPa, Mach number of 2.92, and 
temperature of 92.4K. For all LES/RAS cases, periodic 
boundary conditions were used in the homogeneous- 
spanwise direction. All walls were assumed adiabatic and 
no slip, except for the vertical wall during coarse grid runs, 
where wall functions were applied due to excessively-high 
y+ values observed there during runtime. For all runs, an 
extrapolation condition was applied at the exit plane, and a 
characteristic condition was enforced at the far-field plane. 

The governing equations were integrated using a dual-time- 
stepping approach, in which the diagonalized approximate 
factorization (DAF) method is used for integration in 
pseudo time with a maximum Courant-Friedrichs-Lewy 
(CFL) number of 100.0, and a second-order three-point 
backwards finite difference approximation was used for 
integration in physical time [22,23]. Convergence was 
defined as achieving a residual reduction of 2.5 orders-of- 
magnitude or as reaching a limiting number of 15 sub- 


iterations. The physical time step was set to O.lps for each 
LES/RAS simulation. Following initialization, the flowfield 
was evolved for a minimum of 25 characteristic flow- 
through times before taking statistics. A characteristic flow- 
through time was defined as the time required for a particle 
in the free stream to travel from the comer of the step to the 
point of reattachment of the free shear layer along the 
inclined wall — approximately 0.2ms. Statistical stationarity 
was confirmed by comparing independent sets of statistics 
gathered at intervals of five flow-through times. After 25 
flow-through times, the flowfield was statistically 
stationary, and statistics were gathered for a final five flow- 
through times at a rate of 0.5ps. These statistics were 
subsequently averaged spatially in the homogeneous 
spanwise direction. The RAS turbulence model constants 
were set according to Menter [8], and SGS turbulence 
model constants were set as suggested by Yoshizawa and 
Horiuti [9]. 

5. Results 

The simulations reproduced the major flow phenomena 
observed in the experiment, as seen in Fig. 4; the supersonic 
flow upstream of the cavity detaches at the step corner, 
thereby creating a free shear layer over the cavity. This free 
shear layer attaches along the inclined wall and interacts 
with an oblique shock front standing off the inclined wall. 
Fluid is entrained into the cavity near the location of shear 
layer attachment, and this fluid creates a recirculation zone 
which drives an unsteady low-frequency motion of the free 
shear layer. Note that in Fig. 4, annotations are included 
that indicate major flow features and identify each image’s 
relative location in time, which is given as the value of 
elapsed time, At, divided by the characteristic flow-through 
time, At. In the following subsections, results are presented 
and organized according to study objective, namely: 
assessment of solution accuracy via comparison to 
experiment for the baseline case (R35D), comparison of 
results on the narrow (R8D) and wide (R8DW) grids, 
assessment of grid sensitivity by comparing fine grid 
solutions (R35D and R35L) to those on the coarse grid 
(R8D and R8L), determination of inflow turbulence 
recycling effects by comparing cases with (R8D) and 
without recycling (NR8D), and determination of shock 
sensor effects by comparing cases performed with the 
Ducros switch (R35D and R8D) to those performed with the 
Larsson switch (R35L and R8L). 

Experimental Comparison 

Before comparing results to experimental data, the approach 
boundary layer was confirmed to have developed properly. 
This confirmation was necessary to ensure that the recycling 
procedure did not adversely affect the inflow boundary layer 
statistics. For both the fine grid and coarse grid cases, 
boundary layer velocity profiles at several planes upstream 
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Figure 4 — Instantaneous snapshots of velocity magnitude at the spanwise centerline plane for case R35D with the 
recycle plane, major flow features, and data extraction locations marked. 


of the step were compared to theoretical predictions based 
on Cole’s Law-of-the-Wall formulation and the original 
RAS solutions under the van Driest transformation 
described by Choi, et. al. [10]. These three planes included 
the inflow plane, the donor recycling plane, and a plane at 
the step corner. For each case and at each plane, the 
averaged LES/RAS results agreed well with both the 
theoretical log-layer predictions and with RAS data. The 
comparisons for cases R35D and R8D at the recycling plane 
are included in Fig. 5 to demonstrate the level of agreement. 
In each case and at each plane, the LES/RAS simulations 
reproduced the slope and location of the log-layer well, and 
the blending function anchored the mean LES/RAS 
transition point near the upper edge of the log-layer. Thus, 
the LES/RAS approach and recycling procedure provide 
proper development of the boundary layer upstream of the 
step. 

Mean velocity and static pressure at several planes are 
compared with experimental data in Fig. 6 for each case, 
though details of the baseline case, R35D, are highlighted in 
this section. The results of cases R35L, R8L, and R8DW 


are omitted from Fig. 6 due to their superb agreement with 
cases R35D, R8D, and R8D, respectively, and will be 
discussed later. Note that in Fig. 6, an * denotes a value 
normal to the wall, oo indicates a free stream condition, 8 is 
the approach boundary layer thickness, and all spatial 
dimensions are non-dimensionalized by the cavity depth, D. 
The results of case R35D agree well with the experimental 
data, except within the cavity region. Static pressure along 
the inclined wall and at the chosen planes normal to the 
inclined wall agree well with experiment and mark an 
increase in accuracy, as compared to case NR8RAS. The 
velocity data at planes normal to the inclined wall also mark 
an increase in accuracy. The improved agreement in static 
pressure and velocity data along the inclined wall is no 
surprise, as it is known that RAS models tend to severely 
over-predict the boundary layer recovery rate when free 
shear flows reattach to solid surfaces. However, while case 
R35D more accurately predicts the shear layer spreading 
rate, it is apparent the shear layer is deflected toward the 
cavity floor. Forming a “best-fit” line anchored at the step 
corner that follows the shear layer formed in the R35D case 
results in a deflection angle of 2.85°. 
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Figure 5 — Velocity and blending function in the 
approach boundary layer for cases R35D and R8D at the 
recycling plane. 


Case R8DW was performed to determine whether slow- 
moving, large-scale motions within the cavity were 
artificially limited on the baseline computational domain. 
The results of this case confirmed that no such artificial 
forcing occurred on the original narrow coarse grid. 
Therefore, the authors posit that the experimental data may 
exhibit effects of the tunnel sidewalls and aerodynamic 
fences, which were not modeled in the current study; 
previous efforts have shown that three-dimensional effects 
play a significant role in flow separation even when 
extensive steps are taken in the experiments to ensure some 
level of two-dimensionality [24]. In order to clearly 
illustrate shear layer spreading rate comparisons, each of the 
cavity velocity profiles are shifted vertically in Fig. 7 until 
the point of inflection is approximately aligned with that of 
the experimental data. As expected, fine grid cases predict 
the spreading rate well, and coarse grid cases display 
weaker agreement. Case R35D indicates significant 
improvement as compared to the other cases. It should be 
noted that case NR8D is in less agreement than R8D, which 
indicates an effect of recycled turbulent structure on shear 
layer development and will be discussed later. 


x/D=-1.0 x/D=1.5 


x/D=2.5 x/D=3.5 



Figure 6 — Velocity and pressure for cases R35D (blue dashed line), R8D (magenta dash-dot line), NR8D (orange dash- 
dot-dot line), and NR8RAS (solid red line) to experiment (open triangles) at locations highlighted in Fig. 4. 
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Grid Sensitivity 

Grid sensitivity was assessed by comparing results of cases 
R35D and R8D. Case R8D utilized a once-coarsened 
version of case R35D’s grid in the streamwise and spanwise 
directions; coarsening the transverse direction would have 
yielded unacceptably-large values of y+ for integrate-to-the- 
wall solutions. Based on these two cases, the grid 
sensitivity of the solution is significant. The coarsened grid 
fails to resolve adequate levels of turbulence, as compared 
to the fine grid, which is illustrated in Fig. 8, where 
instantaneous Mach number contours at the spanwise 
center-line are presented. Note that while the turbulent 
fluctuations in the approach flow boundary layer are well- 
resolved for both cases, these fluctuations and 
corresponding turbulent motions are resolved to a lesser 
degree moving downstream of the step on the coarse grid. 


x/D=-1 .0 


x/D=1 .5 x/D=2.5 x/D=3.5 




Figure 7 — Velocity profiles translated vertically to 
demonstrate spreading rate agreement for cases R35D 
(blue dashed line), R8D (magenta dash-dot line), NR8D 
(orange dash-dot-dot line), and NR8RAS (red solid line) 
compared to experiment (open triangles). 



Figure 8 — Instantaneous Mach number snapshots for 
cases R35D (top) and R8D (bottom) at the spanwise 
centerline plane. 


This observation is quantified by examining resolved TKE 
and SGS TKE contributions downstream of the step, which 
are plotted in Fig. 9 at several locations through the shear 
layer. The percentage of resolved turbulence is calculated 
according to Eq. (1), where TKE res and TKE sgs are the 
resolved SGS TKE, respectively. 


%TKE res 


— x 100 

TKE S GS +TKE res 


( 1 ) 


Peak-to-peak values from Fig. 9 are used with Eq. (1). 
Whereas in case R35D, the turbulence is 71%, 80%, and 
85% resolved at each plane, respectively, the turbulence is 
only 30%, 36%, and 50% resolved at the same planes for 
case R8D. Typically, the desired level of resolved 
turbulence for LES of high Reynolds number flows is in the 
80-90% range [25-27]. In case R8D, the dominance of 
TKE sgs leads to more reliance on the SGS model, which in 
this case yields a larger disagreement with the experimental 
data. 



Figure 9 — TKE profiles for cases R35D and R8D. The 
blue dashed line and red solid line are the SGS TKE for 
cases R35D and R8D, respectively, and the gray dash- 
dot-dot line and green dash-dot line are the resolved 
TKE for cases R35D and R8D, respectively. 

Inflow Recycling 

The effect of inflow turbulence recycling on the mean 
flowfield statistics was isolated by rerunning case R8D 
without recycling. The flow upstream of the step was 
simulated using pure RAS, since without the sustained 
coherent turbulent structure within the approach boundary 
layer, the boundary layer was at risk of partially re- 
laminarizing. As can be seen in Figs. 6 and 7, the results of 
case NR8D are significantly different from those of case 
R8D; the shear layer is significantly higher, the spreading 
rate is lower, the static pressure within the cavity is higher, 
and the shear layer reattachment point is further 
downstream. From the pressure data, it is evident that the 
oblique shock makes a shallower angle with the inclined 
wall as compared to case R8D. Explanation of such a 
discrepancy requires consideration of several codependent 
fluid mechanical phenomena, including the transport of 
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mass through the shear layer and the reattachment point’s 
impact on mass entrainment into the cavity. Turbulent 
transport through the shear layer is illustrated in Fig. 10 by 
lineplots of the Reynolds shear stress component, u"v". 
Through the shear layer, significantly less mass transfer 
occurs for case NR8D as compared to cases NR8RAS, R8D, 
and R35D, thereby driving the static pressure within the 
cavity higher. As the shear layer moves vertically, the 
reattachment point moves further downstream, thereby 
reducing the incident angle of the shear layer with the 
inclined wall. As this angle lessens, the amount of fluid 
being entrained into the cavity reduces further. Contrary to 
earlier examinations of this flow, the current results 
demonstrate the criticality of sustaining coherent turbulence 
within the approach boundary layer when employing the 
most recent developments to VULCAN for such a flow 
[5,28]. 



Figure 10 — Reynolds shear stress component for cases 
NR8RAS (red solid line), R35D (blue dashed line), R8D 
(black dash-dot-dot line), and NR8D (magenta dash-dot 
line). 

Shock Sensor 


The effect of shock sensor choice is assessed by comparing 
cases R35D and R35L and cases R8D and R8L. The 
Ducros sensor was considered the baseline for this study. 
Both shock sensors are designed to compare the magnitude 
of the divergence and curl of velocity; however, the Larsson 
sensor includes logic that toggles the interface flux 
approximation to either fully-upwind or fully- symmetric, 
versus the Ducros sensor’s blending of the two. The result 
of this difference is illustrated in Fig. 11, where an 
instantaneous snapshot of the shock sensor is presented for 
cases R8D and R8L. Note how the Larsson sensor tends to 
confine the upwind flux formulation to regions near the 
oblique shock. The Ducros switch uses the centrally- 
differenced flux formulation through most of the cavity 
region and behind the oblique shock. However, the Ducros 
switch employs the upwind flux formulation above the shear 
layer, since by design, any small fluctuation in otherwise 
smooth, inviscid flow triggers such switch behavior. 
Results of cases R35L and R8L are statistically equivalent 


to those of R35D and R8D, respectively. For such reason, 
the results are omitted from Fig. 6, and therefore, the mean 
flow statistics are independent of whether the Ducros or 
Larsson shock sensor is used for the current simulations. 



-2 0 2 x/D 4 6 8 


Figure 11 — Instantaneous shock sensor contours for 
cases R8D (top) and R8L (bottom) at the spanwise 
centerline plane. 

6. Conclusions 

Simulations of a supersonic cavity flow exhibiting a 
reattaching free shear layer utilizing recent developments to 
the VULCAN CFD code were conducted in order to: assess 
the applicability of LES/RAS, isolate the effect of inflow 
turbulence recycling, determine any sensitivity of the 
solution to grid density, and evaluate the effect of shock 
sensor choice for use with a hybridized flux scheme. Two 
grids were used — a fine grid consisting of 35 million cells 
and a coarsened version of the same grid constructed of 8.8 
million cells. A baseline run was performed on the fine grid 
using the Ducros sensor. Two cases were performed on the 
coarsened grid with the Ducros switch toggling the inflow 
turbulence recycling on and off. The cases performed with 
recycling were repeated using the Larsson switch to isolate 
the effect of shock sensor choice. Results utilizing the fine 
grid compared quite well with experimental data outside of 
the cavity region. Through the cavity region, velocity 
profiles indicated excellent shear layer spreading rate 
agreement but disagreement in vertical location of the shear 
layer. This discrepancy may be the result of three- 
dimensional sidewall effects not considered in the current 
study. Even though measures were taken during the 
experiment to promote two-dimensionality, such as 
aerodynamic fences, previous studies have indicated that 
such approaches likely do not completely alleviate the 
effects of three-dimensionality [24]. Based on the 
comparisons of cases R35D and R8D, the solution is 
strongly-sensitive to grid density. Grid convergence was 
not empirically proven based on mean statistics; however, 
simulations performed on the fine grid proved to resolve a 
reasonable fraction of the turbulence, whereas the coarse 
grid simulations exhibited dominance of the SGS modeled 
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turbulence. Cases R8D and NR8D suggest inflow 
turbulence recycling plays a pivotal role in proper shear 
layer development. Without recycling, significantly less 
mass is exchanged through the shear layer, allowing higher 
static pressures to build within the cavity. Subsequently, as 
the shear layer lifts, less fluid is entrained into the cavity, 
and the reattachment point moves further downstream. The 
oblique shock makes a shallower angle with the inclined 
wall, thereby causing severe disagreement in velocity and 
pressure data as compared to experiment. Finally, the 
choice of shock sensor was found to have a negligible effect 
on the mean flow statistics. The sensors of Ducros and 
Larsson were considered, and the results illustrate the 
relative insensitivity to the detailed functional form chosen 
for the shock sensor. Future work will address the three- 
dimensionality concerns raised in the current study. 
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